High-Speed Propulsion of Flexible Nanowire Motors: Theory and Experiments 
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Micro/nano-scale propulsion has attracted considerable recent attention due to its promise for 
biomedical applications such as targeted drug delivery. In this paper, we report on a new experi- 
mental design and theoretical modelling of high-speed fuel-free magnetically-driven propellers which 
exploit the flexibility of nanowires for propulsion. These readily prepared nanomotors display both 
high dimensional propulsion velocities (up to ^ 21/im/s) and dimensionless speeds (in body lengths 
per revolution) when compared with natural microorganisms and other artificial propellers. Their 
propulsion characteristics are studied theoretically using an elastohydrodynamic model which takes 
into account the elasticity of the nanowire and its hydrodynamic interaction with the fluid medium. 
The critical role of flexibility in this mode of propulsion is illustrated by simple physical arguments, 
and is quantitatively investigated with the help of an asymptotic analysis for small-amplitude swim- 
ming. The theoretical predictions are then compared with experimental measurements and we 
obtain good agreement. Finally, we demonstrate the operation of these nanomotors in a real biolog- 
ical environment (human serum), emphasizing the robustness of their propulsion performance and 
their promise for biomedical applications. 



I. INTRODUCTION 

Micro/nano-scale propulsion in fluids is challenging 
due to the absence of the inertial forces exploited by bio- 
logical organisms on macroscopic scales. The diflSculties 
are summarized by Pur cell's famous "scallop theorem" 
[1], which states that a reciprocal motion (a deforma- 
tion with time-reversal symmetry) cannot lead to any 
net propulsion at low Reynolds numbers. The Reynolds 
number, Re = pUL/ ji^ measures the relative importance 
of inertial to viscous forces, where p and p are the den- 
sity and shear viscosity of the fluid, while U and L are 
the characteristic velocity and length scales of the self- 
propelling body. Natural microorganisms inhabit a world 
where Re ~ 10~^ (flagellated bacteria) to 10~^ (sper- 
matozoa) [2, 3], and they achieve their propulsion by 
propagating traveling waves along their flagella (or ro- 
tating them) to break the time-reversibility requirement, 
and hence escape the constraints of the scallop theorem 
[2, 4]. Because of the potential of nano-sized machines 
in future biomedical applications [5], such as targeted 
drug delivery and microsurgery, interdisciplinary eflForts 
by scientists and engineers have recently resulted in ma- 
jor advances in the design and fabrication of artificial 
micro/nano-scale locomotive systems [6-9]. 

Broadly speaking, these micro/nano-propellers can 
be classified into two categories, namely chemically- 
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powered nanomotors [6-8] and externally-powered pro- 
pellers [9]. Chemically-powered nanomotors generally 
deliver higher propulsion speeds, but due to the require- 
ments for chemical fuels and reactions, their applications 
in real biological environments face a number of chal- 
lenges. Externally-powered propellers are often actuated 
by external magnetic fields. Note that these externally- 
powered locomotive systems are often referred to micro- 
or nano-swimmers in the literatures, but strictly speak- 
ing, they do not represent true self-propulsion because 
of the presence of non-zero external torques. In this pa- 
per, we reserve the terminology, "swimmers", to force- 
free and torque-free self-propelling bodies and refer to 
externally-powered locomotive systems as propellers, or 
motors. 

According to their propulsion mechanisms, externally 
powered propellers can be further categorized into three 
groups (see the summary presented in Table I). The 
first group includes helical propellers [13, 14], as inspired 
by helical bacterial flagella [11], which propel upon ro- 
tation imposed by external magnetic flelds. The second 
group of propellers relies on a surface to break the spatial 
symmetry and provide one additional degree of freedom 
to escape the constraints from the scallop theorem, and 
hence are termed surface walkers [15-18]. Finally, the 
third type of propellers, referred to as flexible propellers, 
exploits the deformation of flexible filaments for propul- 
sion. The new nanomotor presented in this paper falls 
into this category. Dreyfus et al. [10] were the first to re- 
alize the idea experimentally by fabricating a 24/im long 
propeller based on a flexible filament, made of paramag- 
netic beads linked by DNA, and attached to a red blood 
cell. Actuation was distributed along the filament by the 
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TABLE I: Comparison between natural micro-organisms and different experimentally realized externally-powered fuel-free 
propellers, divided into three main categories: (a) flexible propellers; (b) rigid helical propellers; (c) surface walkers (i.e. requiring 
a surface for propulsion). We report the maximum dimensional speeds [/max and the maximum dimensionless speeds [/max = 
[//L/, and their corresponding characteristic lengths L and actuation frequencies /. (*: estimated from the dimensionless 
speed in Ref. [10]). 



Type of Propellers 



Schematic/Micrograph 



Maximum dimensional speed 

[/max [M^/s] 



Maximum dimensionless speed 

?7max = [//L/ (XIQ-^) 



Escherichia coli [11] 
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Flexible propeller [10] 
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Surface walker [15, 16] 
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paramagnetic beads; the presence of the red blood ceh 
broke the front-back symmetry, and allowed the propa- 
gation of a traveling wave along the filament, leading to 
propulsion. Recently, Gao et al. [12] proposed a flexible 
nanowire motor made of only metallic nanowires (with 
three segments of Au, Ag and Ni) readily prepared using 
a template electrodeposition approach, and able to swim 
at speeds of up to C/ ~ 6/im/s for a size of L ^ 6.5/im. In 
contrast to the propeller proposed by Dreyfus et al [10], 
the actuation in the device of Gao et al [12] acted only 
on the magnetic Ni portion of the fllament (the head), 
while the rest of nanomotor was passive. 

In the current paper, we present both a new design and 
a theoretical modelling approach for a flexible nanowire 
motor which offers an improved propulsion performance 
(up to U ^ 21/im/s at an actuation frequency (/) of 
35Hz) , approaching thus the speed of natural microscopic 
swimmers, such as Escherichia coli {U ^ 30/im/s at 
/ = lOOHz) [11] while using a lower frequency. The effect 
of size and frequency can be scaled off by nondimension- 
alizing the propulsion speed by the intrinsic velocity scale 
(the product of body length and frequency, Lf) to ob- 
tain a dimensionless propulsion speed, U/Lf^ which can 
be interpreted as the number of body lengths travelled 
per revolution of actuation (or also referred to as the 
stride length in terms of body length in the biomechan- 
ics literature). The nanomotor put forward in this pa- 
per displays remarkable dimensionless propulsion speeds 
compared with natural microorganisms and other artifl- 
cial locomotive systems (see summary of the literature 
and the current results in Table I). 

After presenting the experimental method and its per- 
formance, we study the propulsion characteristics of this 
new high-speed flexible nanomotor theoretically via an 
analytical model. The critical role of flexibility in this 
mode of propulsion is established flrst using simple phys- 
ical arguments, followed by an asymptotic analysis which 
predicts the fllament shape and propulsion speed in dif- 
ferent physical regimes. The theoretical predictions are 
compared with experimental measurements and we ob- 
tain good agreement. The improved propulsion perfor- 
mance of the new fuel-free nanowire motor makes it 
attractive for future biomedical applications, which we 
further illustrate by demonstrating the performance of 
the propulsion mechanism in an untreated human serum 
sample. 



II. HIGH-SPEED PROPULSION 

A. Nanomotor design and fabrication 

The nanowire motors described in this paper were pre- 
pared using a common template-directed electrodeposi- 
tion protocol. In contrast to the previous three-segment 
(Ni-Ag-Au) design by Gao et al. [12], the new design 
relies primarily on a 1.5/im-long Ni head and a 4 /im- 
long flexible Ag tail (see a Scanning Electron Microscopy 



(SEM) image in Fig. lb). A 0.3//m-long Au segment was 
also included (adjacent to the Ni segment) to protect the 
Ni segment from acid etching during the dissolution of 
the Cu sacriflcial layer, and to allow functionalizing the 
motor with different types of biomolecules and cargos. 
Both the Ni and Au segments have a diameter of 200nm. 
While the Ni segment has a length of 1.5/im useful to 
generate sufficient magnetic torques, only a very short 
segment of Au (0.3/im) was used to minimize the over- 
all fluid drag of the nanomotor. Flexibility of the silver 
segment (Fig. lb) was achieved by its partial dissolution 
in hydrogen peroxide solution [12]. The dissolution step 
leads also to hydroxyl products that chemisorb on the 
Ag surface and result in AgOH and Ag2 surface prod- 
ucts. The dissolved Ag fllament had a reduced diameter 
of approximately lOOnm. For the hydrodynamic model 
considered in this paper, the rigid short Au segment is hy- 
dro dynamically indistinguishable from the rigid Ni seg- 
ment, and hence the Ni and Au segments are considered 
in the model as a single rigid 1.8/im-long segment (1.5/im 
Ni-h 0.3/im Au), i.e. the nanomotor has a total length of 
5.8/im. 

The speed of a nanomotor was measured using Meta- 
Morph 7.6 software (Molecular Devices, Sunnyvale, CA), 
capturing movies at a frame rate of 30 frames per sec. 
The trajectory was tracked using a Metamorph track- 
ing module and the results were statistically analyzed 
using Origin software. The speed measured in this man- 
ner is a time- averaged distance travelled per unit time. 
The measurements were performed when the nanomo- 
tors had reached an equilibrium position (in which case 
the image of the nanowire would stay focused under the 
microscope), which leads therefore to the time-averaged 
measurement of U in the laboratory frame. The equi- 
librium distance between the nanomotor and the bottom 
surface was estimated, by varying the focal plane of the 
microscope, to be at the scale of a few microns. The 
detailed experimental procedures can be found in ESIf. 



B. Propulsion performance 

The flexible nanomotors were driven by a magnetic 
fleld with an unsteady component of amplitude iJi, ro- 
tating sinusoidally in a plane perpendicular to a constant 
component, Hq. The magnetic fleld precessed about the 
direction of the constant magnetic fleld at an angular fre- 
quency = 27r/. The nanomotor was observed to propel 
unidirectionally (straight trajectories) in the direction of 
the constant magnetic fleld. In Fig. 2 we show two nearby 
identical nanomotors under the actuation of the external 
magnetic fleld at / = 20Hz (see Video If). These two 
nanowires propel at essentially the same swimming speed 
along the same direction (the red lines are their trajecto- 
ries in a period of 2 seconds), illustrating the stability of 
this mode of propulsion. For helical propellers [13, 14], 
swimming is due to the rotation of rigid chiral objects and 
hence the swimming kinematics scales linearly with the 
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FIG. 1: (a) Schematic representation of a Ni-Ag nanowire motor, and notation for the model, (b) Scanning Electron Microscopy 
(SEM) image showing the topography of Ni-Ag nanowire which was partially dissolved in 5% H2O2 for 1 minute. 




FIG. 2: Two identical nanomotors swimming under the same 
magnetic field at a frequency / = 20Hz. The red lines display 
the superimposed location of the nanomotors over a 2-second 
interval. 



applied field: a reversal of the direction of rotation of the 
magnetic field leads to propulsion in the opposite direc- 
tion for these rigid chiral objects. In contrast, the flexible 
nanowire motors here exhibit uni-directional swimming, 
independent of the rotational direction of the external 
magnetic field. This is due to the nonlinear swimming 
kinematics arising from the flexibility of the nanowire. 
This simple test illustrates the fundamental difference 
between the propulsion of rigid chiral objects and flex- 
ible propellers. In our case, the direction of swimming 
can be controlled by altering the orientation of the axial 
constant component of the magnetic field, Hq. 

We further show in Fig. 5(a) the trajectories of the 
same nanomotor at different frequencies (see captions for 
details) over a 3-second period (see Video 2t). Upon 
the settings Hi = lOG, Hq = 9.5G, and / = 15Hz, we 
are able to achieve a propulsion speed of U = 14.3 ± 
2.46/im/s. The speed of 20 different nanomotors were 
measured, with all other experimental conditions kept 
fixed; the values of the swimming speeds, [/, reported 
in this paper are averaged quantities over these differ- 



ent nanomotors. One meaningful method of comparing 
the propulsion speed between various propeller designs 
consists in scaling the speed with the only intrinsic char- 
acteristic velocity scale of the propeller L/, where L is a 
characteristic body length, and / is a characteristic fre- 
quency. This allows to quantify the distance travelled by 
the propeller in terms of body lengths per revolution of 
rotation. Escherichia coli bacteria [11] typically propel 
with U /Lf ^ 0.03 body lengths per revolution, while the 
flexible nanomotor reported here was able to travel 0.164 
body lengths per revolution at / = 15Hz (see Table I for 
detailed comparison). The maximum dimensional speed 
achieved was U = 20.8 ± 3.07/im/s with / = 35Hz, cor- 
responding in that case to ^ 0.1 body lengths per rev- 
olution. We then experimentally measured the speed- 
frequency characteristics of these flexible nanowire mo- 
tors (results shown as symbols in Fig. 5b). In the next 
section we present a simple physical model for the loco- 
motion of flexible nanomotors, and compare our theoret- 
ical predictions with these experimental measurements. 



III. A MINIMAL MODEL FOR FLEXIBLE 
NANOMOTORS 

A. Chiral propulsion 

In this section, we illustrate the working principles of 
the flexible nanowire motors. First, we establish that it 
is essential for the nanowire to deform in a chiral fashion 
in order to achieve propulsion. 

For low Reynolds number incompressible flows, the 
governing equations are the Stokes equation Vp = /iV^u, 
and the continuity equation V-u = 0, where /i is the shear 
viscosity, and p and u are the fluid pressure and the veloc- 
ity field respectively. Two properties of the Stokes equa- 
tion can be used to deduce the necessity of the nanowire 
being chiral in order to achieve net swimming, as shown 
by Childress [3]. First, it can be shown that the mirror 
image of a Stokes flow is also a Stokes flow. Therefore, 



suppose a nanowire swims with a velocity U along its 
rotation axis, then its mirror image will also swim at the 
same velocity U (see Fig. 3). Second, since time does not 
appear in Stokes equation, it only enters the problem as a 
parameter through the boundary conditions. This leads 
to the time reversibility of the Stokes equation, meaning 
that the velocity field u reverses its sign upon a t ^ — t 
time reversal. In the context of our nanowire motors, 
suppose the nanowire propels at a velocity U, then when 
time is reversed, the nanowire will propel at a velocity 
— U (Fig. 3). If the deformation of the nanowire is not 
chiral, the mirror image of the nanowire can be super- 
imposed with the original nanowire, and the only thing 
reversed in the mirror image is the rotational kinematics 
{i.e. if the original nanowire rotates clock- wisely, its mir- 
ror image will have exactly the same shape but rotates 
counter-wisely; note that the translational velocity is un- 
changed in the mirror image). In this case, one can also 
notice that the kinematics in the mirror image is the same 
as a time reversal of the original kinematics, except that 
the translational velocity is also reversed for the case of 
time-reversal (— U, due to the time-reversibility of Stokes 
flows). In other words, we have now two nanowires (a 
mirror-imaged nanowire and a time-reversed nanowire) 
having exactly the same deformation kinematics but with 
opposite translational velocity (— U = U), and therefore 
we conclude that this can happen only if the translational 
velocity is identically zero (U = 0). Therefore non-chiral 
deformation cannot lead to net propulsion. This simple 
physical argument shows that a combination of rotational 
actuation and nanowire flexibility is critical for this mode 
of propulsion. Recently, the dynamics of tethered elas- 
tic fllaments actuated by precessing magnetic flelds has 
been studied [19-26] and chiral deformation along the 
fllament has been found to produce propulsive force and 
fluid pumping. The swimming behaviours of an unteth- 
ered flexible magnetic fllament displaying chiral deforma- 
tion was also addressed computationally [27]. 
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FIG. 3: Physical explanation of the necessity of chiral de- 
formation in achieving propulsion. If the deformation is not 
chiral, the kinematics of the mirror image of the nanowire is 
identical to the time-reversed kinematics, leading to U = 0. 



(radius a = 50nm, length L = 4/im) as a classical Euler- 
Bernoulli beam [35]. We then employ a local fluid drag 
model, known as resistive force theory [4], to describe 
the fluid-body interaction. The use of a local and lin- 
ear theory signiflcantly simplifles the analysis and is ex- 
pected to provide quantitative agreement because geo- 
metric nonlinearities and nonlocal hydrodynamic effects 
were proven to be subdominant for gentle distortions of 
a slender body in previous work [29, 32, 33]. 



B. Model setup 

Next we show that a simple model taking into account 
the elasticity of the nanowire and its hydrodynamic in- 
teraction with the fluid medium captures the essential 
physics and provides quantitative agreement with exper- 
imental measurements. We flrst solve for the detailed 
shape of the silver fllament, we then predict the propul- 
sion speed, and flnally we compare our results with the 
experimental measurements. Theoretical modelling of 
this type belongs to the general class of elastohydrody- 
namical problems, which has recently received a lot of 
attention in the literature [22, 28-34]. 

Under our theoretical framework, we model the mag- 
netic Ni segment as a rigid slender rod (radius = 
lOOnm, length L^=1.8/im) (the short Au segment is con- 
sidered to be part of the rigid rod in this model, as dis- 
cussed above, see Fig. 1), and the flexible Ag nanowire 



Notation for the model is shown in Fig. 1(a). The 
external magnetic fleld precesses about the z-axis in 
the clock-wise direction, and can be described as H = 
{HicosQt, —HisinQt^ Hq) = i^o(hcosl]t, — hsinl^t, 1), 
where h = Hi/Hq is the dimensionless relative strength 
of the rotating (Hi) and constant (Hq) components of the 
magnetic field. We study the regime where the nanowire 
follows synchronously the precessing magnetic field, ro- 
tating at the same angular frequency (ft) as the magnetic 
field about the z-axis. In addition, we can move in a ro- 
tating frame in which the magnetic field is fixed and the 
shape of the fiexible nanowire does not change with time. 
In this frame, the precessing magnetic field is given by 
H = i^o(h, 0, 1), and the nanowire has a non-changing 
shape r{s) = {r±{s)^z{s)) = {x{s)^y{s)^ z{s)) in a back- 
ground flow, V5, rotating counter-clockwise about the z- 
axis: V5 = l^e^ x = r^(— 7/,x,0), where is the unit 
vector in the z-direction and s is the arclength parameter 
along the fllament. 
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C. Elastohydro dynamics at low Reynolds number 

We describe the fluid-body interaction by resistive 
force theory, which states that the local fluid drag de- 
pends only on the local velocity of the filament relative 
to the background fluid (although in a non-isotropic fash- 
ion). This is thus a local drag model which ignores hy- 
drodynamic interactions between distinct parts of the fil- 
ament, but was shown to be quantitatively correct for 
gentle distortions of the filament shape [29, 32, 33]. The 
viscous force acting on the filament is thus expressed as 

fvis = -K||tt+^^(l-tt)]-U, (1) 

where t(s) is the local tangent vector, u(s) = U — 
is the local velocity of the filament relative to the back- 
ground flow V5, and U is the swimming velocity of the 
nanomotor. Here, and are the tangential and nor- 
mal drag coefficients of a slender filament (L a) and 
are given approximately by 

^11 log(L/a) - 1/2 ' log(L/a) + 1/2 ' ^ ^ 

where fi is the viscosity of the fluid (water, fi = 
10~^Ns/m^). Since the Ni and Ag segments have differ- 
ent aspect ratios {Lm/ctm for the Ni segment), a different 
set of drag coefficients ( ^^y, £,m±) is used for the rigid 
segment. 

When the flexible Ag filament of the nanomotor is de- 
formed, elastic bending forces arise trying to minimize 
the bending energy. This elastic bending force can be 
obtained by taking a variational derivative of the en- 
ergy functional £ = ^ A (^d'^r/ds'^) ds, where A is 
the bending stiffness of the material. The elastic bend- 
ing force is then given by 

felastic — (3) 

Since we are in the low Reynolds number regime, iner- 
tial forces are negligible, and the local viscous fluid forces 
balance the elastic bending forces, fvis+ felastic = 0, which 
yields the equation governing the filament elastohydro dy- 
namics 

[C\\tt + U{l-tt)].u = -A~ (4) 

The flexible Ag filament is clamped to the magnetic 
Ni segment, which is assumed to be rigid and straight. 
Hence, its position vector is given by r^(s) = r \s=l 
+t \s=L {s - L), where s G [L, L -h Lm]- 

D. Nondimensionalization 

We now nondimensionalize the variables and equations 
and identify the relevant dimensionless parameters gov- 
erning the physics of this problem. Specifically, we scale 



lengths by L, rotation rates by 17 = 27r/, times by Q.~^ ^ 
velocities by fluid forces by ^^L^l], fluid torques by 
elastic forces by Ajl?^ and elastic torques by 
AjL. Using the same symbols for simplicity, the dimen- 
sionless elastohydro dynamic equation now reads 

[7-Ht + (l-tt)].u=-Sp-^^, (5) 

where we have defined 7 = ^x/^y, and Sp = 

^(^_L^/^)^^^ is termed the sperm number, which char- 
acterizes the relative influence of the fluid and bending 
forces. 



E. Asymptotic analysis 

The geometrical nonlinearity of Eq. (5) renders the 
elastohydrodynamic equation only solvable via numerical 
simulation in most situations. Here we are able to illus- 
trate the essential physics of flexible nanomotor propul- 
sion analytically via an asymptotic analysis for the case 
where h = i^i/i^o is small. Such an approximation drops 
the geometrical nonlinearities and, as will be shown be- 
low, separates the task of determining the filament shape 
and swimming velocities of the nanomotor, as the axial 
velocities are one order of magnitude smaller than the 
transverse velocities, the axial swimming kinematics be- 
ing thus slaved to the transverse kinematics [34]. Even 
with this simple model, we find that the theoretical pre- 
dictions agree well with the experimental measurements. 
In the experiments, we do not observe very significant 
distortion of the flexible Ag filament, which might ex- 
plain the success of this simple model. 

As the nanomotor was observed to propel unidirec- 
tionally in the z-direction in the experiments (i.e. the 
direction about which the actuating magnetic field pre- 
cesses), we write the swimming speed as U = (0,0, C/) 
and aim at predicting the leading order swimming speed 
in h. We do not expect any O(h^) deformation nor swim- 
ming velocities, and hence the appropriate expansions for 
the deformation of the nanowire and the swimming speed 
are given by 

r^(z) = h r^,(z) + h^r^,(z) + Oi^^), (6) 
t/ = h t/i+h' t/2 + 0(h'), (7) 

where we have 5 ~ z + 0{\i\ 

The elastoydro dynamic equation is a fourth-order par- 
tial differential equation in space, and needs thus to be 
supplied with four boundary conditions. We prescribe 
dynamic boundary conditions at the free end ^ = 0, 
requiring it to be force- free and tor que- free, which is 
d'^Y/dz'^{z = 0) = and d'^Y/dz^{z = 0) = respec- 
tively. Since the deformed shape rotates about the z- 
direction, without loss of generality, we assume the Ni 
head lies on the x — z plane. We then prescribe kine- 
matic boundary conditions at the other end z = 1: 
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FIG. 4: (a) Variation of the dimensionless propulsion speed at second order, U2, with the sperm number, Sp. (b) Superimposed 
snapshots of predicted three-dimensional shape of the Ag nanowire at equal time intervals (t = [T/6,2T/6, ...,5T/6,T] from 
dark to bright color, where T is the period of the rotating magnetic field), for four different sperm numbers. The Ni head is 
not shown here for simplicity. 



r^{z = 1) = (6,0) and dr^/dz{z = 1) = (h,0). From 
experimental observations, the value of b is seen to be 
negligibly small (6^0) and is difficult to measure ac- 
curately. Here, for simplicity, we thus take 6 = in our 
calculations below. In this geometric model, we assume 
that the slope of the magnetic Ni head, dx{z)/dz{z = 1), 
follows the slope of the external field (h = Hi/Hq), which 
is a good approximation when the magnetic field strength 
is strong or when the frequency of actuation is low, such 
that the Ni head can align closely with the external mag- 
netic field. The magnitude of the magnetic torque can be 
compared with the viscous torque acting on the Ni head, 
and their ratio is given by the so-called Mason number. 
Ma. The ratio varies from 0.018 - 0.12, for frequency 
varying form 5Hz to 35Hz. One can also compare the 
magnetic torque to the characteristic viscous torque act- 
ing on the Ag filament, and it varies from 0.13-0.93, for 
the same range of frequency. In both cases. Ma is thus 
typically small and is at most 0(1) at high frequencies. 
Therefore, within the range of frequency explored in the 
experiment, our geometrical model is considered to be a 
valid approximation. At higher frequencies, we would get 
Ma 1, which would play a role in the boundary condi- 
tion at 2: = 1. In that regime, the viscous torque would 
dominate the typical actuation torque by the magnetic 
field, and the Ni segment would therefore not be able to 
align with the magnetic field closely. We expect that the 
slope of the Ni rod might then be smaller than that of 
the magnetic field, and the phase lag between the mo- 
tion of the Ni segment and the magnetic field could be 
substantial. As a result, a degradation in the propulsion 
performance would be expected to occur in this regime. 



1. Determining the flexible filament shape: 0(h) 
calculations 

At order 0(h), the local viscous force is given by 
fvis = h {-yi{z),xi{z),--f~^Ui) + O(h^). From here, 
we can integrate the 0(h) local viscous force in the z- 
direction over the entire nanomotor and since this total 
force needs to vanish because of the absence of exter- 
nal forces, we find Ui =0: swimming occurs therefore 
at order O(h^). The elastic force is given by feiastic = 
-h '&Y>-^{d^xi/dz^, d^yi/dz^, 0) + O(h^). Balancing the 
local viscous and elastic forces in the transverse direc- 
tions yield the hyper-diffusion equations [32] 



-yi Sp" 



xi — Sp 



dz^ ' 
dz^ ' 



(8) 
(9) 



which govern the first order filament shape. The general 
solution to this system of partial differential equations is 
given by 



^l(^) = X] (^P ^riZ) , (10) 

n=l 
8 

^l(^) = X] -^nTn exp (Sp TnZ) , (11) 
n=l 

where is the n-th eight roots of —1, and are com- 
plex constants to be determined by the boundary condi- 
tions. The boundary conditions to this order at 2: = 
are given by d^xi/dz^{z = 0) = d^yi/ dz^{z = 0) = 
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d^xxldz^{z = 0) = d^yxjdz^^z = 0) = 0. The ap- 
propriate boundary conditions ai z = 1 are given by 
xi{z = 1) = yi{z = 1) = 0, dxi/dz{z = 1) = 1, and 
dyi/dz{z = 1) = 0. The 0(h) filament shape is now 
completely determined. 



2. Determining the swimming speed: O(h^) calculations 

At order O(h^), the local viscous fluid force acting on 
the flexible fllament in the z-direction is given by 



- l)A(^) 
l)A(z) - 



< z < 1, 

l<Z<l + i 



where Im = Lm/L, and we have introduced the func- 
tion A{z) = yi{z)dxi/dz{z) — xi{z)dyi/dz{z). Since the 
nanomotor is overall force- free, the second order swim- 
ming speed U2 can be determined by integrating the local 
viscous fluid in the 2:-direction over the entire nanomotor 
and requiring this total force to vanish, i.e. 



I 



' fvis2^^ — 



(12) 



and we see that the swimming speed is slaved to the 
flrst order fllament shape {xi{z),yi{z)) via the function 
A(z). Upon simpliflcation with Eqs. (8) and (9) and the 
boundary conditions at 2; = 0, we obtain 



U2 = 



•7 



dxi d^xi 



(13) 



dz dz^ 



d^xi 

dz^ 



Oyi O^yi 
dz dz^ 



dz^ 



where a = Crn\\/^\\- dimensional form, the leading 
order swimming speed, is given by 



U 



mil J 



(14) 



dxi d^xi 
dz dz^ 

+ 0(h'). 



d^xi 

dz'' 



+ 



d'yi 



dz dz^ 



d'yi 

dz'' 



As in previous work [22, 23, 33, 34, 36, 37], we observe 
that this mode of propulsion relies on the drag anisotropy 
of slender fllaments, 7 = ^±/^|| ^1. Indeed, when 7 = 1, 
the swimming speed vanishes. Note that for very slender 
fllaments, we have 7 2 (see Eq. 2). 

We also see that the swimming speed scales quadrati- 
cally with the relative strength of the rotating and con- 



stant components of the magnetic fleld, U 



for 



h ^ 1. This scaling is conflrmed by a complementary 
asymptotic calculation valid for low sperm numbers (see 
Appendix for details). 

Next, we plot our predicted dimensionless second or- 
der swimming speed as a function of the sperm num- 
ber Sp (Fig. 4a) together with the predicted fllament 



shapes (Fig. 4b), and observe three different character- 
istic regimes. The sperm number Sp is the most impor- 
tant dimensionless group governing the propulsion per- 
formance. For Sp <C 1, bending forces dominate and the 
fllament is effectively straight (Fig. 4b-i). Hence, the fll- 
ament motion is almost kinematically reversible, and it 
produces small propulsion. Quantitatively, a small Sp 
asymptotic analysis presented in the Appendix reveals 
that the dimensionless swimming speed grows with the 
fourth power of the sperm number, C/2 ^ Sp^, for Sp <C 1. 
On the other hand, from Eqs. (10) and (11), we see that 
most deflection is concentrated around a small region 
< z < 1/Sp, when Sp ^ 1, due to the exponential de- 
cay of the solution amplitude. In this regime (Sp ^ 1), 
the viscous forces dominate, and propulsion is ineflicient 
because a large portion of the fllament has small de- 
flection and thus experiences drag but contributes to no 
thrust (Figs. 4b-iii & iv). As a result, we expect optimal 
swimming to occur when Sp is of order one, where the 
total drag of the nanomotor is kept low while the drag- 
induced bending is fully exploited to produce propulsion. 
This is conflrmed in our calculation, and we observe the 
optimal sperm number to occur at Sp ~ 2, which gives 
a maximum propulsion speed of U2 ~ 0.042 (Fig. 4a). 
The fllament shape close to optimal swimming (Sp = 2) 
is shown in Fig. 4b-ii. 



F. Comparison with experimental measurements 

Under flxed magnitude of the rotating and constant 
components of the magnetic fleld, the swimming speed 
of a nanomotor was measured with the frequency of the 
magnetic fleld varying between to 35Hz. The exper- 
iment was repeated on the same nanomotor for three 
different settings of magnetic fleld strengths (shown us- 
ing three different symbols with error bars in Fig. 5b; a 
total of 20 different nanowires were sampled). The ro- 
tating magnetic fleld strength Hi was kept constant at 
Hi = lOG, and the constant magnetic fleld strength was 
set to be Hq = 14. 3G (blue squares), Hq = 11. 8G (red 
circles), and Hq = 9.5G (green diamonds). 

We then compare in Fig. 5b our theoretical predictions 
(solid lines) with experimental measurements, plotted as 
swimming velocity vs. frequency (main flgure) or sperm 
number (inset). In our theoretical model, the value of the 
bending stiffness A of the flexible fllament is unknown. 
Standard bending stiffness of pure silver is inapplicable 
here since the dissolution of silver in hydrogen perox- 
ide rendered the fllament a porous structure with signifl- 
cantly reduced strength and a different chemical compo- 
sition (Ag20, AgOH). A value of A = 3.6 x lO-^'^Nm"^ 
flts, with the least total squared errors, the experimental 
data with Hi = 10G,i7o = 14. 3G (blue squares), which 
is the case where our model is expected to work best as 
the ratio h. = Hi/Hq is the smallest. This bending stiff- 
ness is then used to predict the propulsion speed under 
different magnetic fleld settings (green and red solid lines 
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FIG. 5: Dependence of the nanomotor swimming speed on the actuation frequency, (a) Superimposed trajectories of the same 
Ni-Ag nanomotor at different frequencies / = 5, 10, 15, 30Hz (as indicated) over a 3-second period (red hues), with Hi = lOG 
and Ho = 9.5G. The scale bar is lO/xm. (b) Speed-frequency characteristics of flexible nanowire motors. Symbols represent 
experimental data for different setups of the magnetic fields: blue squares {Hi = lOG, Ho = 9.5G); red circles {Hi = lOG, Ho 
= 11.8G); green diamonds {Hi = lOG, Ho = 14. 3G). Error bars show standard deviations of the measured speeds (20 samples). 



The solids lines show the theoretical predictions (Eq. 14) with A 
of the swimming speed on the sperm number, Sp. 



3.6 X 10 Nm . The inset in (b) displays the dependence 



in Fig. 5b, see captions for details). 

The theoretical model is seen to capture both qualita- 
tively and quantitatively the speed-frequency character- 
istics of these flexible nanomotors. Qualitatively, the rate 
of change of the swimming speed with respect to the fre- 
quency increases at low frequencies (/7 ~ for small /, 
as shown in the Appendix) , but then gradually decreases 
as the frequency continues to increase, and eventually 
levels off at high frequencies. Physically, increasing the 
actuation frequency is equivalent to increasing the sperm 
number. When the frequency is varied from to 35Hz, it 
corresponds to a variation of the sperm number Sp from 
to 2.6, experiencing a degradation in swimming perfor- 
mance beyond the optimal sperm number Sp ~ 2, which 
corresponds to a frequency of around 15Hz in our exper- 
iment. This degradation manifests as a less-than linear 
speed-frequency variation (since the dimensional swim- 
ming speed scales as LQ^ linearly in Q) beyond the fre- 
quency 15Hz, resulting in the level-off at higher frequen- 
cies. As noted above, at very high frequencies, the mag- 
netic Ni head will be unable to follow synchronously the 
rapid rotating magnetic field. The dynamics of propul- 
sion will be more complicated in that regime, and the 
simple model presented here will likely be inapplicable. 

The agreement between our theoretical model and our 
experimental results is very satisfactory. The discrepan- 
cies are larger for the setup Hi = lOG, Hq = 9. 50 (green 
lines and squares), which is expected because h ^ 1.1 
in this case and the asymptotic assumption of small h is 



less valid. Note that our measurements did not sample 
the low Sp regime as in our experiments, swimming at 
low frequencies appear to be significantly influenced by 
Brownian motion. 

Our model has only one fitting parameter, the bend- 
ing stiffness A, which - as explained above - we fit to the 
bottom data set in Fig. 5b, and use to predict the other 
two data sets. The estimated value we obtain from the 
fitting is equivalent to a pure silver nanowire of diame- 
ter ^ 6nm (with elastic modulus, E = 80OPa), which 
is much smaller than the diameter of flexible segment 
observed. This is expected because the chemical compo- 
sition of silver is altered after the dissolution, and a large 
portion of the flexible nanowire is indeed a thick layer of 
surface byproducts formed after the chemical reaction, 
which contributes little, if any, to the bending strength. 
The diameter of the actual structural filament that bears 
the bending loads is difficult to measure experimentally 
(see details of the structure in Fig. lb). In addition, non- 
uniform chemical reactions lead to strong local defects or 
points of weakness along the nanowire, which might sig- 
nificantly reduce the bending strength. We can compare 
our estimated bending stiffness, A = 3.6 x 10~^^Nm^, 
with the bending stiffness of typical flagella of natural 
microscopic swimmers, such as eukaryotic spermatozoa, 
which also rely on the flexibility of flagella for propulsion. 
These biological filaments have their bending stiffnesses 
ranging from IQ-^^Nm^ [38] to IQ-^^Nm^ [39], which is 
the range in which our estimated value lies. 



10 



















f=Os 


f=0.5s 


\ 

f = 1s 


f= 1.5s 



FIG. 6: Time lapse images (time as indicated) of the motion 
of nanowire motor (velocity, U = 15/im/s) in human serum 
at / = 15Hz, with Hi = lOG and Ho = 9.5G. Scale bar is 
5/im. 



IV. DISCUSSION 

In this work, we designed and fabricated a high-speed 
fuel-free nanomotor utilizing the flexibility of nanowires 
for propulsion. These flexible nanomotors demonstrate 
a number of advantages: first, the fabrication process 
is relatively simple and involves a common template- 
directed electrodeposition protocol of nanowires; second, 
these nanowire motors are able to propel at high speeds, 
both dimensional (up to ~ 21/im/s) and dimensionless 
(up to 0.164 body lengths per revolution), and their per- 
formance compares very well with natural microorgan- 
isms and other synthetic locomotive systems; third, they 
are actuated by external magnetic field and do not re- 
quire specific chemical environments (fuels) for propul- 
sion, which is preferable for biomedical applications. In- 
deed, the performance of the nanomotors reported here 
is not affected by the presence of ions or other chem- 
ical species, and they are able to propel equally well 
in real biological settings. As an illustration, we have 
placed these flexible nanomotors in human serum, and 
observed similar propulsion behaviours (Fig. 6 and Video 
Sf). This demonstrates an exciting potential of these flex- 
ible nanomotors for future biomedical applications such 
as targeted drug delivery systems, or cell manipulation. 

The fundamental physics of the flexible nanomotors 
has been illustrated by a simple analytical elastohydro- 
dynamic model. The propulsion characteristics were ex- 
perimentally studied and compared with the theoretical 
model, with good agreement. Strictly speaking, the re- 
sults of the asymptotic model presented in this paper are 
valid only for h <C 1. However, as shown in other pre- 
vious studies which compared asymptotic results with 
numerical and experimental studies [29, 32, 33], these 
asymptotic models often remain valid even up to h 1, 
meaning that geometrical nonlinearities do not play very 
significant roles. Our results also ignore the hydrody- 
namic effect of the bottom surface close to which the 
nanomotors are propelling. As we estimated experimen- 
tally the distance of the filaments to the surface to be on 
the order of microns, and therefore on the order of the 
swimmer size, we do not expect strong hydrodynamic ef- 



fects from the surface, which might explain the success 
of our simple modelling approach. Further progress in 
theoretical modelling most likely have to be obtained nu- 
merically. More accurate yet complicated descriptions 
of the hydrodynamic interactions can be achieved using 
methods such as slender body theory [40] , or regularized 
Stokeslets [41]. Tension, self-spinning, and twist strains 
of the filament may also be considered for improvements. 
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Appendix A: Swimming at low sperm numbers 

In this appendix, we consider another physically in- 
teresting asymptotic limit, the low sperm number limit, 
F = Sp^ <C 1. The results in this asymptotic limit are 
not expected to provide quantitative agreement with the 
experimental measurements, since the value of F in the 
experiment is typically large. Nevertheless, this analy- 
sis still allows us to reveal different scaling behaviours 
of the propulsion speed at low sperm numbers. Keaveny 
and Maxey [27] investigated the propulsion of a flexi- 
ble filament with distributed magnetic actuation. With 
a resistive force model, they considered the low sperm 
number limit and found that the propulsion speed scales 
with the fourth power of the sperm number, U ~ Sp"^, at 
low sperm numbers. We will follow closely their method 
of solution and perform similar calculations here to show 
that the same scaling holds for our flexible nanowire mo- 
tor subject to boundary actuation. An explicit formula 
for the leading order propulsion speed in F will be de- 
rived. Expanding this formula for small h will confirm 
our scaling of U with h from the small-h asymptotic anal- 
ysis in the main text. 

First, the problem is formulated below, taking into ac- 
count the effects of twisting, self-spinning, and inexten- 
sibility of the nanowire. Denoting N(5, t) and M(s, t) 
as the resultant internal force and moment on a cross 
section, the local force and moment balances are 

^ = (^ii-a)(t-u)t + ^^u, (Ai) 

^ + t X N = ^ra'^{iv • t)t, (A2) 

where uj{s^t) is the angular velocity, and = 47r/i is the 
resistive coefficient for the viscous torque produced by 
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self-spinning (rotation about its own local axis, t) of the 
filament. The internal moment M(s, t) has a constitutive 
relation 



M = Atx ^+K,^t, 

OS OS 



(A3) 



where Kf is the twist modulus of the filament and ^(s, t) 
is the twist angle. In contrast to the propeller studied by 
Keaveny and Maxey [27], the magnetic torque does not 
come into the local moment balance in our case, but only 
through the boundary condition. The boundary con- 
ditions are given by the balance of external forces and 
torques at the ends of the flexible filament. We have a 
free end at s = 0. The external forces and torques at 
s = 1 are given by the total viscous force and viscous 
torque together with the magnetic torque acting on the 
Ni segment, which is modelled as a slender rigid rod: 



N{s = 0) 
M{s = 0) 
N(5 = 1) 



0, 
0, 



-Lm [(Cm|| - ^m_L)(t • u)t + ^ 



M{s = 1) = Mt 



=1 xH 



ml. 



{^txu+^[a;-(t.a.)t] 



(A4) 
(A5) 

(A6) 



S = l 

(A7) 



where M = M^a^TrL^ is the strength of the magnetic 
moment of the Ni segment and Mg = 485 x lO^A/m is 
the spontaneous magnetization of Ni. 

To study the low sperm number limit, we adopt the fol- 
lowing nondimensionalizations: we scale times with 
lengths with L, H with Hq, elastic forces with A/L'^, and 
elastic torques with A/L. With these nondimensionaliza- 
tions, the dimensionless equations (using the same vari- 
ables for simplicity) now read 



dH 



l)(t •u)t + u] 



RT((jj •t)t = tx— -+txN 

OS"^ 



^ ds \ ds^ 



where R = irO^/i^L'^, K = Kt/A, and V 
dimensionless boundary conditions are 



N(5 
M(5 



0) 
0) 



0, 
0, 



N(. = 1) = -r/3 x I [(7-1 - l)(t . u)t + u] + l^u; X t| 



M(s = 1) = C^t 1,^1 xH-r^ 



{ 



(A8) 
(A9) 
Sp''. The 

(AlO) 
(All) 

2 

(A12) 



2 3 



characterizing the relative strength of the magnetic and 
elastic torques. Finally, we have the inextensibility con- 
dition 



t(s,i) •t(s,i) = 1. 



(A14) 



Eqs. (A8) through (A14) completely describe the full 
swimming problem (within the realm of resistive force 
theory and classical elastic beam theory) without mak- 
ing any assumption. There is no restriction on the va- 
lidity of the solution to this system, but the solution 
has to be obtained numerically, with special attention on 
the nonlinearities arising in the differential equation and 
the boundary conditions. To make analytical progresses, 
we consider the asymptotic limit F ^ 1, and calculate 
the leading order swimming speed in V. Following the 
method and notations by Keaveny and Maxey [27], we 
assume the filament attains a constant shape at steady- 
state and rotates about the 2;-axis synchronously with the 
external magnetic field, hence we write the dimensionless 
steady-state conformation of the filament as 

x{s, t) = -b{s) cos[t + (/)(5)], (A15) 
y{s, t) = -h{s) sin[t + (/)(s)], (A16) 
z{s,t) = a{s)^Ut, (A17) 



where U = U/LQ is the dimensionless swimming speed 
in the z-direction, whereas a{s)^ K^)^ ^.nd (j){s) are ge- 
ometrical unknowns to be determined. Like the small-h 
asymptotic analysis in the main text and in Ref. [27], 
here we have considered unidirectional swimming in the 
2;-direction. In addition, since the functions a{s)^b{s) 
and (l){s) are independent of time, we only find the solu- 
tion at one specific time, t = [27]. The magnetic field 
at t = is given by H = (h, 0, 1). We seek expansions in 
F in the form of 



NiF + 0(F2), 
MiF + 0(f2), 
da I 

ds 

dbi 



;-(t. 
(A13) 



;)t] 



N = No- 
M = Mo 

da dao 

ds ds 
db dbo 

ds ds 
d(j) d(j)i 

ds ds 
dip dip- 

ds ds 



ds 

f o(r2), 
-r + o(r2), 



r + c>(r2), 
^ + o(r2), 



(AI8) 

(A19) 
(A20) 

(A21) 

(A22) 

(A23) 
(A24) 



where Im = L^/L, j3 = ^m±/^±, Im = Cm±/'?m||, and 
Cm = MsO^TrLmHoL/A is a dimensionless parameter 



and similar expansions hold for other variables. With 
these expansions, we can express the local tangent t(s) 
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and velocity u(s, t) as 

t(5) = to + tir + o(r2), 

ds ^ ^ ds 



with boundary conditions 



(A25) 
(A26) 



+ l-^,-^,^)r + o(n 
m{s) = (o, -60, u^) + (^001, t/i) r + o{v^). 

(A27) 

In the fonowing section, we wiU perform the calculations 
order by order. 



1. 0(r°) calculations 

The O(r^) local balance of forces and torques are given 



by 



ds 



0, 



(A28) 



= to X No + to X ^ + K— ( ^to ) , (A29) 
with boundary conditions 



No(s 


= 0) 


= 0, 


(A30) 


Mo(s 


= 0) 


= 0, 


(A31) 


No(s 


= 1) 


= 0, 


(A32) 


Mo(s 


= 1) 


= Cmto \s=l XH, 


(ASS) 



and inextensibility condition 

The solution at this order is given by 

No{s) = 0, 

dao 1 

ds ~ + h^ ' 

dbo — h 

a/i + h^ ' 

*o(s) = i/o = 0. 

2. 0(r) calculations 



(A34) 

(ASS) 
(AS6) 

(AS7) 
(ASS) 



The 0(r) local balance of forces and torques are given 
(7-^-l)(to-uo)to + uo, (AS9) 



by 
dNi 



ds 



= to X Ni + to X — ^ + ti X — ^ + K—^to, 

(A40) 



Ni(s = 0) = 0^ 
Mi(s = 0) = 0, 
Ni(s = 1) = - 



(A41) 
(A42) 



(7„^ - l)(to • uo)to + uo + Y^o X to 



s=l 

(A4S) 



Mi(s = l) = C™ti U=i xH 

_ r.,2 J to X Uo I ^,1 
P'-n 



[u>o - (to • Wo)to] 



and the inextensibility condition 

dai j^^^i 
ds ds 



(A44) 



(A45) 



From the solution at 0(r°), Eq. (AS7), we can inte- 
grate to get h(j{s) = — hs/x/r+l? + Ci, where C\ is a 
constant to be determined. By satisfying the equations 
and boundary conditions at this order, we find that 



h[l + /3?™(2 + ^„)] 



2Vm?(l+/3/m) 
The solution at this order is then given by 
d{ho<j)i) hs^ Cis^ 



ds 24\/l + h^ 6 
dbi dai d'^i 
ds ds ds 



Bi, 



0, 



(A46) 

(A47) 
(A48) 



where 
Bi = 



24C™ (1 + h^) (1 + pU) 



{ 



(A49) 



by 



2 + CmVl + i^[l + Plm{3 + 2lm)] 
+ 2plm [4 + lm{Q + «m(4 + Plm))]^- 

3. C>(r2) calculations 

The 0(r2) local balance of forces and torques are given 



<m_ 

ds 



2 1 

(7~ - l)(ti • Uo + to • ui)to + ui, 

(A50) 



R{u:i • to + a;o • ti)to = to x Ng + tj x Ni (A51) 

as^ as^ ds'^ 
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with boundary conditions 



N2(s = 0) = 0, 
M2(s = 0) = 0, 



(A52) 
(A53) 



where 

All = A22 = 



h{j-' - 1) , h/3Z^(7-i - 1) 



l + h^ 



l + h^ 



, (A64) 



N2(s = l) = -/ 



(7^^ - l)(to • Ui + ti • uo)to 



+ Ui + (u?o X ti + U?i X to) 



(A54) 



M2(s = 1) = C^t2 |.=1 X 



In 



H-^Z^|^(to 



Ui + ti • Uo) 



A12 = - 

A21 = - 
$1 



1 + h^ 



1 + h^ 



+ [ct^i - (ti • cjo + to • ct;i)to - (to • a;o)ti] 

(A55) 



From the local balance of force (Eq. A50) with the 
boundary conditions at 5 = (Eq. A52), we find that 
N2 = (A^2x, N2y, N2z) is given by 



\/r+i? 



1 + h^ 



^(1) + 



— + 1+ ^ 



1 + h^ 



(A65) 
(A66) 



N2.{s) = 



Vl + h^ 



(7-^ - 1) 



X UlS 



yiTl? Vl + h 
+ H{s) + C3S, 

N2y{s) = —blS, 



h X h ^ 

H(s) + ^ C3S 



vT+i? 



(A56) 
(A57) 



N2z{s) : 



J(s) + + 



Vl + h^ Vl + h 



h X h ^ 



^1+1? 



(A58) 



where Cs is an unknown integration constant and we 
define the functions 



F{s) 



ds 



■Si, (A59) 



24Vr+l? ^ 

G(.) = / F{s')ds' = - ^ + 

^ 120V 1 + h^ 

(A60) 



^0 



720Vl + h^ ^ ^ ' 

(A61) 



^^(5) = / bo{s')F{s')ds\ 
Jo 



(A62) 



Examining the force components N2X and N2z^ we have 
two unknowns, namely Ui and C3. These two unknowns 
are determined by applying the boundary conditions at 
s = 1 (Eq. A54), yielding a 2 x 2 system of equations 



'All Ai2\ (Ui 
iA2i A22) [Cs 



h$i + ^lmG{l) + ^F(l) + H{1) 
$1 

(A63) 



yr+1? 



Vl + h' 



yr+ir 



=^(1) 



:G(l) + 6o(l)F(l) 



(A67) 



Upon solving this system of linear equation, we arrive 
at an explicit formula for the leading order swimming 
speed 



(A68) 



1440 (1 + h^)'^' (1 + Plm^i^plm + 7m) 

|577m[l + /3/m(4 + 3/^)]^ - Ajplm [4 + /3/^(13 + 9/^)] 
- 57m - JmPlm [24 + 30/^ + /3lm (28 + 84/^ + 45/^)] |r 



o(r2). 



(A69) 



Again, we can verify that when we have isotropic drag 
7 = 7m = 1, then no swimming is possible, Ui = 0. 



4. Variation with the relative magnetic field 
strength, h 

First, one can see that Cm, which is the ratio of the 
characteristic magnetic torque to the characteristic elas- 
tic torque, does not enter the formula for Ui (Eq. A69), 
meaning that the absolute value of the magnetic field 
strength or the dipole moment strength has not yet 
played a role in the swimming speed at low sperm num- 
bers. However, the relative strength of the rotating and 
constant magnetic field, h = Hi/Hq, has an interesting 
effect here. From the small-h asymptotic analysis in the 
main text, we knew swimming occurs at O(h^) (Eq. 14), 
and hence the swimming speed (both dimensional or di- 
mensionless) scales quadratically with h, U ~ h^, for 
h ^ 1. This is confirmed by examining Eq. (A69), where 
we have U - hV(l + h^)^/^ = h^+0(h^) when expanded 
for small h. When the dimensionless swimming speeds 
from the two asymptotic analyses are plotted against h 
for small sperm numbers (Fig. 7, different lines repre- 
sent results at various sperm numbers, see the caption 
for details), we see an excellent agreement when h is 
small, illustrating that the swimming speed does increase 
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FIG. 7: Variation of the dimensionless swimming speed, 
U/LQ, with the relative magnetic field strength, h, for differ- 
ent sperm numbers based on the low-Sp calculations: Sp = 0.2 
(dark blue solid line), Sp = 0.3 (light green solid line). The 
red dotted (Sp = 0.2) and black dash-dotted (Sp = 0.3) lines 
are the corresponding results from the small-h calculations. 
Inset: Same plot for Sp = 0.5 (dark orange solid line) and 
Sp = 0.6 (light blue sohd line). The green dotted (Sp = 0.5) 
and brown dash-dotted (Sp = 0.6) lines are the corresponding 
results from the small-h calculations. 



quadratically with h for small h (the dotted lines are the 
small-h asymptotic results). From the low sperm number 
results (solid lines), the swimming speed then experiences 
a maximum when h continues to increase (the small-h re- 
sults are no longer valid in this regime), and eventually 
decreases with further increase in h. From the analytical 
expression (Eq. A69), we see that the maximum swim- 
ming speed occurs at h = ^/2 ^ 1.41. 



5. Variation with the sperm number, Sp 

From the small-h asymptotic analysis in the main text, 
we have illustrated the dependence of the dimensionless 
swimming speed on the sperm number, Sp (Fig. 4). Here, 
via the low sperm number asymptotic results (Eq. A69), 
we see quantitatively that the dimensionless swimming 
speed scales as the fourth power of the sperm number, 
U = U/m - r - Sp^, for low Sp. Since Sp"^ ex /, 
this also means that the dimensional swimming speed 
scales quadratically with the frequency ~ , for small 
/. We confirm this result by plotting the variation of 
the dimensionless swimming speed with the sperm num- 
ber (Fig. 8), and the variation of the dimensional swim- 



ming speed with the frequency (Fig. 8 inset). We com- 
pare the low-Sp asymptotic results (h = 0.2, dark blue 
solid line; h = 0.4, light green solid line) with the corre- 
sponding small-h asymptotic results h = 0.2 (red dotted 
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FIG. 8: Variation of the dimensionless swimming speed, 
U/LQ, with the sperm number, Sp, based on the low-Sp cal- 
culations: h = 0.2 (dark blue solid line), h = 0.4 (light green 
solid line). The red dotted (h = 0.2) and black dash-dotted 
(h = 0.4) lines are the corresponding results from the small-h 
calculations. Inset: Variation of the dimensional swimming 
speed U with the frequency / based on the low-Sp calcula- 
tions: h = 0.2 (dark blue solid line), h — 0.4 (light green 
solid line). The red dotted (h = 0.2) and black dash-dotted 
(h = 0.4) lines are the corresponding results from the small-h 
calculations. A bending stiffness of A = 3.6 x 10~^^Nm^ is 
used in the speed-frequency plot (inset). 



line) and h = 0.4 (black dash-dotted line), so that the 
small-h asymptotic assumption (h <C 1) is expected to 
be valid. We see the results from the two asymptotic 
analyses agree with each other for sufficiently low sperm 
numbers (Fig. 8), or frequencies (Fig. 8 inset), confirming 
the scaling U/LVt - Sp"^ ([/ - Z^), at low sperm numbers 
(low frequencies). 

To summarize, in this appendix, with the help of a 
low-Sp asymptotic analysis, we have confirmed the scal- 
ing U for small h in the main text, and established 
the complementary scalings U/LQ ~ Sp^ for small Sp, 
or equivalent ly U ^ f'^ for small /. Note that the results 
in this appendix are valid for very small Sp numbers, 
and hence are not expected to be useful for experimen- 
tal comparison. In addition. Cm is assumed to be 0(1) 
throughout the calculations. Since Cm cx Im, the results 
here are invalid for <C 1 and /^n ^ 1- 
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